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We discuss the scale-free property of Wilson's numerical renormalization group(NRG) for the 
Kondo impurity problem. The single-particle state of the effective Hamiltonian with a cutoff A is 
described by the wavepacket basis having the scale free property; The energy scale of the system can 
be controlled by the lattice translation of the wavepacket basis with no reference of rescaling of the 
lattice space. We also analyze the role of the Kondo interaction in the context of wavepacket basis 
and then discuss the scaling and renormalization of the Kondo coupling. In addition, we clarify the 
role of the edge state in the lowest energy scale of Wilson NRG. 
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I. INTRODUCTION 

Renormalization group (RG) has been one of the most 
fundamental concept in statistical physics to explain hi- 
erarchical structure of the energy scale in nature^ In 
particular, such a real-space RG approach as block- 
spin transformation provides an intuitive view for coarse 
graining of the fluctuation field and rescaling of the 
length unit; It is well-known that the block-spin trans- 
formation for the triangular lattice Ising model is of- 
ten discussed as a textbook example 2 . For the low- 
dimensional quantum RG, however, the situation is dif- 
ferent; The conventional block-spin transformation of- 
ten fails in extracting the low-energy physics of the one- 
dimensional(lD) quantum system, and improvement of 
the real-space RG has been an attractive frontier in the 
quantum physics. In modern view, this is mainly because 
the block-spin transformation breaks the quantum entan- 
glement associated with the large quantum fluctuation^ 

Historically, two real-space RG approaches not based 
on the block-spin transformation made great success in 
the ID quantum physics: Wilson's numerical renormal- 
ization group (NRG) 5-7 for the Kondo impurity problem 8 
and density matrix renormalization group(DMRG)^i£. 
An interesting point is that these two methods have com- 
pletely different characters. Wilson NRG can directly 
deal with the metallic system with a magnetic impu- 
rity, through mapping to the effective ID chain with a 
boundary. On the other hand, we should remark that 
DMRG includes no scale transformation and no con- 
cept of the running coupling constant. This implies that 
DMRG should be considered as a variational optimiza- 
tion method for the matrix product wavefunctionii - — , 
rather than a RG based on the coarse graining and the 
scale transformation. Thus, Wilson's pioneering work 
has been the peculiar RG method that succeeds in such 
energy-scale control of the massless excitation. What is 
its key mechanism? Answer to this question provides 
an essential insight for the quantum RG. Here, we just 
note that, recently, the entanglement renormalization in- 
terestingly reproduces the correct critical index for ID 
quantum systems^. 



In this paper, we will show that Wilson NRG has the 
scale free property, which is the key of systematical se- 
lectivity of the arbitrary energy scale around the Fermi 
surface. For this purpose, we directly focus on the re- 
lation between the wavefunction and the cutoff param- 
eter A due to the log-discritization. In actual NRG it- 
erations, a free fermion site is added step by step with 
decreasing the energy scale by A and the effective ID 
tight-binding model having the exponential modulation 
of the hopping parameter is eventually generated. The 
qualitative role of this cutoff can be viewed as a regula- 
tor for the infrared divergence, which introduces effective 
length scale £ ~ 1/lnAji 7 - In actual numerical computa- 
tions, however, the cutoff value is usually adopted to be 
A = 1.5 ~ 2.5, which is beyond perturbation level to the 
original Kondo problem(A = 1). The success of Wilson 
NRG with the nonperturbative value of A suggests that 
the Wilson's cutoff scheme should involve some plausible 
mechanism to maintain the nature at A = 1; We clarify 
that the one-particle eigen- wavefunction is modified from 
the plane wave into the wavepacket by the cutoff, and its 
lattice translation explains the selectivity of the energy 
scale. We also find that there is an edge state in the low- 
est energy scale near the Fermi surface. On the basis of 
the wavepacket basis above, we perform the exact diago- 
nalization of the Wilson effective Hamiltonian including 
the Kondo impurity. We then discuss that the scale fac- 
tor of the groundstatc wavefunction should be renormal- 
ized from the value of the non-interacting wavepacket by 
the Kondo interaction. 



This paper is organized as follows. In section II, we 
define the effective tight-binding model and then explain 
its scale free property. In section III, nature of the 
wavepacket is analyzed. In section IV, we resolve the role 
of the Kondo iteration in the context of the wavepacket 
basis. Section V is devoted to the summary and discus- 
sion. 
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II. MODEL AND SCALE FREE PROPERTY 

In Wilson NRG, the essential point is to consider the 
ID lattice fermion model mapped from the degenerating 
free electrons around the Fermi surface. We thus start 
with the Wilson-type Hamiltonian of the spinless free 
fcrmions with the exponentially modulated hopping 

jv-i 

KX= E eA "( C Ul C "+4cn + l), (1) 
n=l 

where c n is a fermion annihilation operator at nth site 
and N denotes the number of sites. We have also intro- 
duced A = exp(A) > 1 for later convenience. Thus n = 1 
corresponds to the smallest energy scale and n = N does 
to the impurity site with the largest energy scaled. Al- 
though, in the original work , the hopping parameter has 
a supplemental coefficient and n in e Xn term takes a half 
integer, the essential physics is the same as Eq. ([1]). 

Let us write the one-particle state as = 
En i J ( n ) c n\fy ■ Then one-particle Schrodinger equation 
in the bulk region is 

e-^(n-l) + f(n|l)=&^(n). (2) 

Note that Eq. t^fy is invariant under the transformation, 
ijj(n) —> (—l) n ijj(n) and E — > — E, which clearly rep- 
resents the particle-hole symmetry. Thus, we basically 
consider the positive energy solution. 

Since the system has no explicit translational sym- 
metry, we employ numerical diagonalization of Eq. ([3]) 
rather than the usual Fourier analysis, for a finite but 
sufficiently large system. We assume the free boundary 
condition and thus what we deal with is the tridiagonal 
matrix. Figure 1 represents the absolute value of the 
one-particle spectrum for A = 0.1 and N — 200, where j 
indicates the label of the eigenvalue in increasing order. 
The Fermi surface is located between j = 100 and 101, 
and j < 100 represent negative energy eigenvalues. Thus, 
the parity in Fig. 1 corresponds to the particle-hole sym- 
metry. We also show the amplitude of the wavefunctions 
corresponding to j = 101, 130 and 160th eigenvalues, in 
Fig. 2. 

The most important behavior in Fig. 1 is that, as 
was already mentioned in Ref£, E basically exhibits the 
exponential dependence E cx ±exp(Aj). We call this re- 
gion of the exponential dependence as "bulk", since the 
corresponding wavefunctions are localized in the bulk re- 
gion of the chain, as can be seen for j = 130 and 160 in 
Fig. 2. On the other hand, we can see that some eigenval- 
ues near the Fermi surface j ~ 100 deviate from the bulk 
lines. We call these states as "edge states", since they 
correspond to the edge modes near the Fermi surface, as 
is the wavefunction of j =101 in Fig. 2. 

We analyze the bulk part of the spectrum in detail. 
In connection with the exponential dependence of the 
eigenvalues, an essential information can be found in the 
bulk wavefunctions for j = 130 and 160 in Fig. 2. The 
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FIG. 1: (Color online) One-particle eigenvalue spectrum of Eq. 
(J2j) . The horizontal axis j indicates the label of the eigenvalue 
in increasing order. 
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FIG. 2: (Color online) One-particle wavefunction tp(n) of j = 
101, 130 and 160. The horizontal axis means the cite index n. 
The bulk states 4>(n) of j — 130 and 160 can be overlapped 
with each other by a lattice translation, while ip(n) of j = 101 
is the edge state. 

primary notable point is that these wavefunctions have 
a very similar wavepacket-like shape, in contrast to the 
usual plane wave for the uniform chain. The localization 
of the wavefunction can be qualitatively understood as 
follows. If a particle carrying a certain energy goes in 
the larger n region, the particle can not excite the larger 
energy bonds, into which the wavefunction can not pen- 
etrate. While the particle goes to the smaller n region, 
the bond of the smaller coupling can not carry the total 
energy of the particle and thus the wavefunction decays 
very rapidly. 

Another important property of the bulk wavefunctions 
is that they can be overlapped with each other by the 
lattice translation; In Fig. 2, the wavefunctions of j =130 
and 160 have the very similar shape. Indeed, we can 
verify that overlap integral of the two wavefunctions after 
the lattice translation is unity within the computational 
accuracy. In order to see this property in analytic level, 
we introduce 

^(n) = e- Xn t 2 (f>{n). (3) 
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and then rewrite Eq. ([2]) as 

<f)(n - 1) + <j>{n + 1) = Ee- Xn cj)(n), 



(4) 



where E = Ee x l 2 . Now let us consider the effect of the 
lattice translation. By shifting n = nl + to and rewriting 
cj){n' + to) — > 4>{n'), we obtain 



1) 



1 



Ee 



(5) 



which becomes identical to Eq. (j4]) by redefining E' = 
Ee~ Xm . In other word, the lattice translation of the 
wavefunction has one-to-one correspondence to the eigen- 
state of the e _Am -scaled energy eigenvalue, which has 
been mentioned in RefJ^ in the operator level. Thus, 
once a certain energy eigenvalue is obtained, the other 
eigenvalues and the corresponding eigenvectors can be 
obtained by the lattice translation combined with E' = 
e~ Xm E. In addition, we note that, if A is sufficiently 
small, we can always find a certain m such that E 1 is 
unity, except for E being exactly zero. This implies that 
the Wilson Hamiltonian is energy- scale- free theory and 
thus the problem can be physically well-defined, although 
the Hamiltonian includes very wide range of the energy 
scale, A ~ A N . 



III. WAVEPACKET BASIS 

We analyze features of the bulk wavepacket itself. For 
this purpose, we concentrate on the positive energy solu- 
tion of Eq. (g]) and consider its "continuum limit" X 0. 
Writing x = An, we have 



\ 2 (j)"(x) + (2- Ee~ x )(t)(x) =0, 



(6) 



where we have omitted C(A 3 ). Note that the translation 
x — > x+a combined with Ee~ a — > E keeps Eq. (j6j) invari- 
ant. Thus we can set E = 1 in the following analysis. It 
should be also remarked that the continuum limit A <C 1 
can be taken with no reference to the lattice space. 

Introducing y = exp(— x/2), we can convert Eq. (J6j) to 
the modified Bessel type equation, 



■ ■ ld_ j_ -_ . , 

dy 2+ ydy + A 2 ^ 2 ' 



2 
V 



4>{y) = o. 



(7) 



Assuming that the boundary is at infinity with respect 
to x, and taking account of the normalizability of the 
wavefunction in i <C — l(y ^> 1), we obtain 4>(y) 
cons: x Ki2 C (cy) with c = \/2/A, 20 where Ki U {z) is the 
modified Bessel function of the imaginary order v. Using 
the formula J °° dzzK-i V (z)Ki V (z) = 2 ^ , we obtain 
the eigenfunction of Eq. ^ as 



i>(x) = f^^ce-^K t2c (ce-^). (8) 

The behavior of this solution can be resolved in special 
limits. When a; > l(i/ < 1), we can easily see that Eq. 
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FIG. 3: (Color online)Amplitude of the wavefunction at the 
impurity site ipj(N) of TV = 50 for A = 0.05, 0.1, 0.2 and 0.5, 
where j indicates the label of an eigenstate. We can see the 
clear exponential dependence on j for the bulk wavefunction. 
The small hump around Fermi surface(j ~ 25) corresponds to 
the edge states of the lowest energy scale. The hump region 
of the edge states extends, as A becomes small. 



(j6]) reduces to the uniform tight binding model, which 
yields 



(x) ~ cxp(— icx). 



(9) 



This is consistent with the asymptotic behavior of the 
modified Bessel function K iu (z) ~ 2 w ~ 1 T{iv)z~ ll/ in the 
z — > limit. Thus, the wavefunction ip(x) exponentially 
decays in x S> 1 with the oscillation. This behavior is 
consistent with the numerical result in Fig. 2, where we 
have checked that the tail of the wavepacket decays as 
e~ A "/ 2 . 

In the limit z 3> 1, the asymptotic form is given by 
Kiv(z) ~ y^exp(— z). For x -C 0(y ^> 1), thus, we can 
see 



exp(- 



ce 



-x/2 



), 



(10) 



which decays very rapidly in x < 0. This is also in good 
agreement with the numerical result in Fig. 2. From 
these, we can see that the localization of the wavepacket 
is basically good and the width of the packet is ex- 
actly described by the decay rate of the high-energy tail, 
£ = 2/A(= 2/ In A), because of Eq. ©. Also this good 
localization of the packet justifies that the effect of the 
boundary can be negligible in the bulk region. 

Since the Kondo impurity is located at the right edge 
of the chain, the amplitude at n — N is of particular 
importance. Here we return to the lattice problem and 
write the index of the wavepacket as j . Figure [3] shows 
the absolute value of ipj (N) , for which we can see the ex- 
ponential dependence of ipj (N) with respect to the index 
j. Because of the resolution of the figure, we plotted the 
case of N = 50. Note that the particle-hole symmetry 
of ([1]) gives a., = — otN-j+i with j < N/2(N =even), 
which can be also confirmed in Fig. [3J According to Eq. 
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©, the tail of the single fermion wavepacket is described 
by the damping factor e~ A ™/ 2 , when the packet is suffi- 
ciently away from the impurity. Thus, we can write the 
amplitude of the wavefunction at the impurity site as 



a 3 exp(-Aj) 
ctj exp(A(j — N 



for 1 <j< f , 
1)) for f + 1 < j < N, 



(11) 

where oij is a coefficient of order of unity for j-th 
wavepacket. Note that otj is a constant in the bulk 
region. However, we should also remark that, as A 
becomes small, the small hump region extends around 
Fermi surface^' ~ 25); The states in the hump region 
correspond to the edge states of the lowest energy scale, 
which will be discussed in the next section. 



IV. EDGE STATE 

In this section, we would like to discuss the bound- 
ary effects in Wilson NRG, which can be seen around 
the Fermi surface(i? = 0) corresponding to the left edge 
of the chain. Figure |4] shows the magnification of the 
eigenvalue spectrum around E = 0(j ~ 100), where the 
scale of the vertical axis is linear. As was mentioned, the 
energy eigenvalue distribution in the bulk part has the 
exponential dependence, which is represented as broken 
curves. However, we can see that the spectrum in the 
vicinity the Fermi surface deviates from the exponential 
dependence and becomes rather linear with respect to 
index j. These states have the lowest energy scale near 
the Fermi surface, which implies that the property of the 
edge state is important. Indeed we can see that the wave- 
function of j = 101 in Fig. 2 localizes at the edge region. 
Although the shape of this wavefunction itself is similar 
to the right side of the bulk wavepacket, the eigenvalue 
distribution exhibits different structure from the bulk re- 
gion. 

In order to analyze the edge states, we consider another 
continuum limit, which may be similar to the usual ap- 
proach for the fermion with the linear dispersion. We 
rewrite <j>(n) — e ±47rn / 2 ry(n), where the phase factor ba- 
sically corresponds to the right or left going Fermi wave 
numbers of the uniform case. Then, equation (j4)) becomes 

T i[r)(n - 1) - r}{n + 1)] = Ee~ Xn T]{n). (12) 



Writing x = An again and maintaining the leading term 
of 0(A), we have T Ee~ x )r](x) = 0. The solution 

of this equation can be easily obtained as 



T}(x) = Cexp(=F*— e - *), 



(13) 



where C is a certain constant. For x 3> 1, r\ — ?>const and 
thus V ~ e ~ x / 2 e ±l7TX / 2X . For x -C — 1, however, r\ be- 
comes highly oscillating, and if its period becomes smaller 
than the lattice space, this continuum limit breaks down. 
Moreover, the wavefunction |-0(a;)| 2 ~ e~ x — > oo in 



x — > — oo. In this sense, the wavefunction (| 1 3|) is pro- 
hibited in the bulk limit. If the system has the left edge, 
however, such a breakdown can be avoided. When the 
boundary condition is tp(0) = ij}{N + 1) = 0, the energy 
eigenvalue is determined as 



E = 2ir\e- x/2 q, 



(14) 



where q — ±1/2, ±3/2 • • •— is a quantum number mea- 
sured from the Fermi surface. In fact, we have E = 
0.2988 for A = 0.1 and q — 1/2, which is in good agree- 
ment with the numerical result 0.2989 in FigHJ As E 
increases, the oscillation of Eq. (|X3|) becomes significant 
and then the wavefunction crosses over to the bulk wave- 
function ©, away from the edge. Since the phase dif- 
ference between n — and 1 in Eq. (jT3j) is about E/2, 
the criterion of the crossover is estimated as E ~ E ~ 2, 
which is also consistent with Fig. @] 

Although the lowest energy behavior is described by 
the boundary excitation (TT4"|) , the role of the edge state 
attracts less interest so far. One reason for this is that the 
practical value of the cutoff is A 2 ~ 1.5 ~ 2.5, for which 
the level space becomes too sparse to detect the linear 
region in the spectrum. Here, it should be noted that 
the generalized Wilson NRG scheme with much smaller 
A(= 1.02) exactly extracts the low-energy spectrum of 
the critical XXZ chain*^ Another reason is that, as in 
Fig. [U the edge wavefunction smoothly varies its shape 
into the bulk wavepacket and the tail of the edge wave- 
function shows the same dumping |i/>| 2 ~ e~ x as the bulk 
wavefunction. Thus, one can increase number of sites 
seamlessly in NRG iterations, with no explicit reference 
of the edge state. Finally, we make a comment on the 
right edge corresponding to the highest energy scale; In 
contrast to the left edge above, the universal eigenvalue 
structure can not be seen for the right-edge states. 




FIG. 4: (Color online)The linear scale plot of the eigenvalue 
spectrum around the Fermi surface for A = 0.1 and N = 200. 
Solid symbols denote the jih eigenvalue. The Fermi surface 
corresponds to E — 0, which is indicated by the horizontal 
dotted line. Solid curves represent exponential dependence of 
the eigenvalues in the bulk region and the broken line indi- 
cates the edge mode (fTi")) . 
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V. EFFECT OF IMPURITY SPIN 

On the basis of the nature of the wavepacket, let us 
discuss the effect of the Kondo impurity in Wilson NRG. 
Assume an impurity spin S with the Kondo coupling Jk 
at the right edge of the chain(n = N). The free electron 
Hamiltonian 7i\ is diagonalized in the wavepacket ba- 
sis, H\ = Ejfjtrfjer, where fj a is a fermion operator 
corresponding to the j-ih eigenvalue of (HJ, and the spin 
index a is noted explicitly. In the following, the eigen- 
state described by fj t „ is called as "orbital electron" for 
convenience. As was discussed in Sec. II, Ej oc exp(2Aj) 
for 1 < j < N/2, and E 3 oc exp(A(iV - 2jj) for 
N/2 + 1 < j < N 22 . On the other hand, the matrix 
element associated with the impurity is determined by 
the tail of the wavepacket at the right edge, 

n inip = j K e XN 

(15) 

where a is the 5 = 1/2 spin(not Pauli) matrix associ- 
ated with the orbital electron, 5* represents the impurity 
spin of S = 1/2, and a takes er Gf,4- ■ Here, it should 
be recalled that, in Wilson NRG, the energy scale of the 
Kondo coupling itself is scaled up as J^e AAr , in the pro- 
cess of the iterative increase of the free electron sites. 

It is instructive to write the matrix elements of "Hi mp in 
the orbital bases. Using Eq. (fTTj) . we have the diagonal 
parts with respect to the spin as, 

(i,a,S\'Hi m p\j,a,S) = ±-^-aiajC itj , (16) 

where 

, e X(N-i-j) > l<i,j<K 

c J e^" 1 ), l<i<f, f + l<j<N 
° w " | e^-V, f + 1 < i < N, 1 < j < f [ 
{ e A(-iV+ i +,--2) ! N + 1 <iJ< N 

Also, the off-diagonal part is given by 

(htA\n imp \j,i,i\) = (i,lA\n imp \jAA) 

= -f"'". <"',;• (is) 

where "f|",JJ- represent the impurity spin. Here we should 
remark that the damping factor of the wavepackets near 
the Fermi-surface ~ N/2) is canceled out with the 
scale factor of the Kondo coupling. This implies that the 
fermionic 2fcF-oscillation of Eq. ([§]) is correctly recov- 
ered. Thus, the phase shift of the Fermi-level state due 
to the impurity scattering can be correctly taken into 
account, although the wavepackets are resolved with re- 
spect to the energy scale. 

The eigenvalues of the impurity term Himp are just 
given by — ?>e XN Jk/^ for the singlet, e XN Jk/^ for triplet 
and the zero energy with the huge degeneracy. By the 
perturbation of Ti.\, the most of the degenerating zero 



eigenvalus are lifted to have a "dispersion" , whose maxi- 
mum energy scale is of order of Em ~ e XN . For Jk > 1, 
we can naively treat H\ perturbation against to 

Himp, since "band width" of the free orbital electrons is 
smaller than the singlet-triplet gap of Himp. In the con- 
text of physics, an orbital electron localizes at the edge 
and forms a singlet or triplet pairs with the impurity spin. 
This is the local Fermi liquid fixed point. 

For the weak coupling case(J/f <C 1), the dominant 
part of the Hamiltonian is described by 1-L\, and %i mp 
is a perturbation for small N . As is well known, how- 
ever, the low-energy physics should flow to the strong 
coupling limit for sufficiently large N. For simplic- 
ity, we assume 1 < i,j < For 2Af-site lattice 
translation 2 - 2 , then, the scale factor Cij in (fl6| and 
([TBI is scaled by Ci+M,j+M exp(2AM), while a^j in the 
bulk region is independent of Thus we can read 
J K e XN ^{N)^(N) = J K e XN ^* +M (N)^ 1+M (N) with 
Jk = Jk exp(2AM). The role of Kondo coupling of a 
Jk at a certain energy scale e Aj can be viewed as that 
of the scaled Kondo coupling Jk exp(2AM) in a higher 
energy scale e x ^ +M \ In the following, we call this fac- 
tor exp(2AM) as "canonical scale factor" of the Kondo 
coupling. 

The above scaling argument might suggest that the 
low-energy physics basically flows to the strong coupling 
limit by the canonical scale factor e 2XM . However, the 
actual spectral flow is rather complicated; Since H\ has 
the same scale factor as Hi mp as well, the competi- 
tion of "Ha and Himp yields the nontrivial Kondo en- 
ergy scale. Of course, it is well established that the en- 
tanglement between the impurity and orbital electrons 
rapidly develops around the Kondo energy scale given by 
exp(—l/2Jx )/\J2Jk- In the present context, this sug- 
gests that the site dependence of the effective Kondo cou- 
pling should be deviated from the naive canonical scale 
factor. 

In order to see the effect of the Kondo interaction, we 
perform exact diagonalization of the Wilson's effective 
Hamiltonian up to 14 orbitals with the subspace of elec- 
tron number N and total-.? 2 = l/2(including the impu- 
rity spin). We use A = In 2/2 ~ 0.347(A = y/2) in order 
to see clearly the role of the impurity within the exact 
diagonalization. Then, we focus on the nature of the 
groundstate wavefunction, rather than the spectral flow. 
Let us write the f-electron density of the groundstate for 
iV-sites system as = (fj+fj-t), where j indicates the 
orbital index. Figure EJa) shows the Jff-dependence of 
Pff ■ Note that, for N — 12, the Fermi surface of the 
free electron is located between j = 6 and 7. As Jk in- 
creases, the f-electron density of j = 7 approaches 0.5, 
which implies that the impurity spin is screened and then 
the remaining spin- 1/2 chiefly moves to the orbital near 
the Fermi surface. At the same time, we can see that the 
density for j ^ 7 orbitals is also modified from the case 
of free electrons, even though they are embedded deeply 
away from the Fermi surface. This implies that the high 
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FIG. 5: (Color online) (a) Kondo coupling dependences of f- 
electron density for AT = 12 with A = In 2/2 ~ 0.347(A = \/2), 
where j denotes the index of the orbitals. The Fermi surface 
of the free orbital electron corresponds to j — 6.5. (b) Size 
dependences of ^-electron density for A ~ 0.347(A = \/2). 
The curves indicate the ^-electron density at the Fermi surface 
Pn/2+11 f° r N = 6, 8, 10, and 12 from down to up. The solid 
curves of represent the f-electron density at the highest 
energy scale orbital. 

energy orbitals are certainly entangled with the impurity 
spin. 

In Fig. Hp}), we show the size dependence of the f- 
electron density at j — N/2 + 1 for N = 6, 8, 10 and 

12. We can easily see that, as Jk increases, , 2 + 

approaches p ~ 0.5; the S z — 1/2 at the impurity is 
passed on the orbitals at the lowest energy scale around 
the Fermi surface, and the number of the orbital sites 
is effectively reduced to be N — 1. Moreover, as N in- 
creases with Jk fixed, p^^i/ 2 f approaches the value of 
the strong coupling limit. This behavior is of course 
equivalent to the iterative process of the Wilson NRG. 
On the other hand, the f-electron density in the highest 
energy scale p^) shows little size dependence, which is 
almost negligible within the scale of FigfSJb). 

In Fig. |51 we show the real-space distribution of the 
spin density of the groundstate for N = 14. Note that 
the charge distribution of all Jk is uniform in the real 
space, indicating the spin-charge separation. The singlet 
pair of the impurity and the orbital electron itself is in- 
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FIG. 6: (Color online)Spin distribution in the real-space rep- 
resentation for = 14 and A = In 2/2. The open(solid) circle 
indicates Jk = 3.0(0.5). n denotes the lattice index, where 
n = 14 corresponds to the impurity site. The dotted line 
means the spin distribution of the free orbital of j = 8. 



visible to a spin expectation value. However, one can 
extract an essential information from the distribution of 
magnetization passed on the orbital electrons from the 
impurity site. We first discuss the strong coupling case 
Jk = 3.0. We have confirmed that the spin distribution 
for Jk — 3.0 is basically identical to that of Jk = 100. 
The dominant single-particle component in the ground- 
state is the orbitals around the Fermi surface(j = 7 and 8 
for N — 14). On the other hand, the groundstate wave- 
function also includes certain amplitudes of the higher 
energy orbitals, as can be seen in Figs. [5j The decay rate 
of the spin density in the strong coupling limit is expected 
to be close to the canonical value A = In 2/2 ~ 3.47, but 
to have a certain deviation from it. In Fig. H2 the open 
circle at n =even is very similar to the dotted line of the 
free orbital of j = 8, but a precise analysis for the decay 
rate of the spin density for Jk — 3.0 at n =even sites 
actually gives A* = 0.312 ~ 0.9A. In addition, we can 
see that the spin density at ri=odd sites clearly deviates 
from the single particle result. The above results of the 
S z distribution illustrate that the wavefunction in the 
strong-coupling limit includes the nontrivial correlation 
effect, where A is slightly renormalized to A* ~ 0.9A. This 
suggests that the scale factor could be also renormalized 
to e 2A M from the canonical value in the wavefunction 
level. 

We next turn to the case of Jk = 0.5, which is in the 
crossover regime of the Kondo coupling within N = 14. 
For Jk = 0.5, the total-iS z of the orbital electron is 0.233, 
which implies that the screening rate of the impurity spin 
is 47%. Then, we can see that the decay rate reduces to 
be A* = 0.215 ~ 0.62A, which illustrates that the wave- 
function in the crossover regime is entangled with the 
states in the wide range of the energy scale. In the con- 
text of the wavepacket basis, this implies that hybridiza- 
tion of the higher-energy orbitals generates an eigenvalue 
of the lower-energy scale. As the system size increases, 
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we have confirmed that the renormalized A* approaches 
the strong coupling limit value 0.312. In this sense, the 
size dependence of A* might be related with the scaling 
function of the Wilson NRG(Eqs. VIII.45, 46 or IX.58 in 
Rcf. 5 ). For the detailed analysis of A*, the wavefunction 
of much larger system size is needed beyond the exact 
diagonalization. 

VI. SUMMARY AND DISCUSSION 

We have clarified that the key mechanism of Wilson 
NRG is the scale free property of the Wilson Hamilto- 
nian JT]); The cutoff A results in the wavepacket basis 
described by Eq. ([5]) and its lattice translation enables 
for us to control the energy scale of the system with no 
reference to rescaling of the length unit. In addition, we 
have shown that the Hamiltonian (JXJ) has the edge states 
at the lowest-energy scale, which has not been mentioned 
so far. On the basis of wavepacket, we have performed 
the exact diagonalization including the Kondo coupling. 
We then revealed that the nontrivial effect of the Kondo 
coupling can be detectable as a renormalization of the 
decay rate of the groundstate spin density; The value of 
A* is close to the canonical value A in the strong coupling 
limit, while the renormalization effect becomes large in 
the crossover regime. Since the effective length scale of 
the ground state wavefunction is given by £* = 1/A*, it 
may be efficient to increase the number of the basis in 
the crossover regime. 

In the context of the uniform ID critical system, 



the critical phenomena are usually characterized by the 
power-low decay of physical quantities. However, the ex- 
ponential modulation introduces the cutoff of the infrared 
divergence and then the criticality of the uniform system 
would be converted to the nontrivial renormalization of 
A*. In particular, the connection of the renormalized 
scale factor in Wilson NRG to the DMRG analysis of 
the Kondo screening cloud for the uniform chain with a 
boundary- 2 ^ is an interesting problem. We also remark 
that, as in Refj 1 ^, the regulator for the infrared diver- 
gence due to such scale-free modulation of the interaction 
is a generally applicable idea to the ID quantum many 
body system. For further research, it may be essential to 
clarify the direct relation between the entanglement and 
the renormalized scale factor. Recently, the connection 
of Wilosn NRG to the matrix product state is pointed 
out in Refs . 24 ' 25 . Moreover, very recently, the bound- 
ary critical phenomena can be correctly captured by the 
entanglement renormalization! 26 ' 27 We believe that the 
present result stimulates a new frontier in the quantum 
RG. 
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